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The one-dimensional (ID) isotropic frustrated ferromagnetic spin-1/2 model is considered. Classi- 
cal and quantum effects of adding a Dzyaloshinskii-Moriya (DM) interaction on the ground state of 
the system is studied using the analytical cluster method and numerical Lanczos technique. Cluster 
method results, show that the classical ground state magnetic phase diagram consists of only one 
single phase: "crural". The quantum corrections are determined by means of the Lanczos method 
and a rich quantum phase diagram including the gapless Luttinger liquid, the gapped chiral and 
dimer orders is obtained. Moreover, next nearest neighbors will be entangled by increasing DM in- 
teraction and for open chains, end-spins are entangled which shows the long distance entanglement 
(LDE) feature that can be controlled by DM interaction. 



in 

(N 



CO 



I 

C 

O 

o 



> 

(N 

in 
o 



PACS numbers: 75.10.Jm; 75.10.Pq 



I. INTRODUCTION 

An important goal in study of quantum spin systems is 
the search for novel phases emerging from competing in- 
teractions between particles. Low dimensional quantum 
spin systems typically exhibit strongly correlated effects 
which drive them toward new regimes with no classical 
analog. Many properties of the systems in these regimes 
or "quantum phases" can be understood if we explore 
their ground state and low-lying excitations. 

Therefore a significant fraction of current research on 
such systems has focussed on understanding of frus- 
trated systems, which exhibits a variety of exotic quan- 
tum states^. In particular, the frustrated spin chains are 
described by competing interactions between the nearest- 
neighbor (Ji) and next-nearest-neighbor (J2) interac- 
tions. The Hamiltonian of the model is given by 



iY 



H = 



.7=1 



(1) 



where Sj represents the S — 1/2 operator at j-th site of 
the chain. The model with both antiferromagnetic inter- 
actions Ji , J2 > is well studied^— . Also the model 
(fTj) with ferromagnetic and antiferromagnetic interac- 
tions ( J\ < 0, J2 > 0) (frustrated ferromagnetic model) 
has been a subject of many studies^— . It is known 
that the ground state is ferromagnetic for a = tj^t < j. 



At a c 



the ferromagnetic state is degenerate with a 



singlet state. The wave function of this singlet state is 
exactly know n 15 ' 16 . For a > \ , the ground state is an in- 
commensurate singlet. It has been long believed that at 
a > i the model is gaples o 7 ' 1 ^ but the one-loop renormal- 
ization group analysis indicate d 13 ' 19 that the gap opens 
to a Lorentz symmetry breaking perturbation. However, 
existence of the energy gap has not been yet verified 
numerically.^ 3 -. Using field theory considerations it has 



been proposed^ 7 - that a very tiny but finite gap exists 
which can be hardly observed by numerical techniques. 

In vicinity of a — j, the singlet ground state energy 
behaves as Eq = (a — jY , where /3 is a critical exponent. 
Using variational approaches^, it has been shown that 
the quantum fluctuations definitely change the classical 
critical exponent 7 -^ and yields j3 = | which has been 
confirmed numerically 2 ^. 

Beside a general interest in understanding frustrations 
and phase transitions in the model systems described 
by the Hamiltonian in Eq(l), it helps people to under- 
stand intriguing magnetic properties of a novel class of 
edge-sharing copper oxides, described by the F-AF frus- 
trated model 2 ^— . Recently some novel magnetic prop- 
erties were discovered in a variety of quasi-one dimen- 
sional materials that are known to belong to the class 
of Dzyaloshinskii-Moriya (DM) (D.(Si x Sj)) magnet to 
explain helical magnetic structures. The relevance of an- 
tisymmetric superexchange interactions in spin Hamilto- 
nians which describe quantum antiferromagnetic systems 
was introduced phenomenologically by Dzyaloshinskii^ 5 .. 
Moriya showed later, that such interactions arise nat- 
urally in the perturbation theory due to the spin-orbit 
coupling in magnetic systems with low symmetry 2 ^. 

Some multifcrroics cuprates, such as LiCuVO ^ 27 ' 28 , 
LiCu-iO^ and Cu-iGeO^ are expected to be described 
by DM interaction. This has stimulated extensive inves- 
tigations of various properties which are created by the 
DM interaction. However, it is difficult to handle the 
DM interaction analytically and interpret experimental 
data. A Numerical analysis then helps us to understand 
experimental observation and even expand our knowl- 
edge about many interesting quantum phenomena of low- 
dimensional quantum magnets. In the present work, we 
address a different problem in the subject of the frus- 
trated ferromagnetic spin-1/2 chains. We consider the ID 
spin-1/2 frustrated ferromagnetic model with added DM 
interaction and study the classical and quantum mag- 
netic ground state phase diagram of the system. 
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The outline of the paper is as follows. In the next 
section the classical cluster method will be outlined and 
phases will be obtained. In section III we present our nu- 
merical results of the exact diagonalization calculations 
on the ground state properties of the model. In section 
IV, the entanglement between different spins will be in- 
vestigated by calculating the concurrence function. Fi- 
nally we conclude and summarize our results in section 
V. 





II. CLASSICAL PHASE DIAGRAM 

The Hamiltonian of the ID spin- 1/2 frustrated ferro- 
magnetic model in presence of a uniform DM interaction 
is defined as 



N 

0=1 



J 2 Sj.Sj +2 



D.(SjXS j+1 ] 



(2) 

where D = Dz. This Hamiltonian shows a model with 
fully broken spin rotational symmetry. It is useful to be- 
gin with simple classical considerations that shed light on 
the possible ground state of the fully quantum mechani- 
cal problem. When S is large, it is adequate to substitute 
Sj = S(sin9j cos 4>j, sin#j sin (j)j, cosOj) in Eq. ^ and 
find solutions that satisfy dH/d8j = and dH/d(f)j = 
for all j. But, here we use the rather unknown cluster 
metho d 30 i 31 of Lyons and Kaplan (LK method). Briefly 
recall that method with assuming periodic boundary con- 
dition. Then we can easily see that Eq. @ can be rewrit- 
ten as: 



H c — 2J h c (Sj-i, Sj , Sj+i 



(3) 



where the "cluster energy" involve three neighboring 
spins is 

h c {S l7 S 2 ,S 3 ) = ^Ji(Si.S 2 + S 2 .S 3 ) 

+ D.(Si xS 2 + S 2 xS 3 )} 



J 2 Si.Ss 



It is clear that 



h c ( Sj-i, Sj, 5j+] 



(4) 



(5) 



so easily one can find the minimum of h c . If the cor- 
responding state propagates, i.e. if there is a state of 
the whole system such that every set of three successive 
spins gives the minimum h c , then according to Eq. (|5|), 
this state will be a ground state of the Hamiltonian. This 
is the main idea of the LK cluster method as applied to 
the present problem. This method has also been applied 



FIG. 1: (color online). The schematic picture of angles made 
by the spins in a cluster of three spins. 



for models with higher dimensions or with open bound- 
ary conditions^. 

Now let's minimize h c . Because the DM vector is se- 
lected along the z axis, we consider only coplanar (spin 
dimensionality d=2, i.e. XY spins) state. Without loss 
of generality label the angles made by the end spins with 
the central spin 9, 9' (FigfTJ). So the cluster energy is 
became 

h c (6,e') = -cos0-cos0'+7(sin0-sm0') 
8 ^ 

+ 2acos((9 - 9')}, (6) 

where 7 = jj^j. Differentiating gives the conditions for 
stationarity 

\Ji 



J-^i j sin 6 + 7 cos 9 -2a sin(0 - 9') j 
Hp j sin 6' - 7 cos 9' + 2a sin(0 - 6') } 



= 
= 0. 



(7) 



Before going through the Eq. ^ , we check two simplified 
cases as follow. 

First case: We have put D = and the frustrated fer- 
romagnetic Heisenberg model recovered which well stud- 
ied numerically and analytically. In this case, solutions 
are: 



{9,9') 
{9,9') 

cos 9q = — 



(0, 0), (tt, 7r)(0, n){n, 0) {I sing - type) 
{9q, —9q) {chiral — type), where 
1 

4o 



(8) 



The {tt, tt) solution, which leads to the ordinary antifcr- 
romagnetic state, is never lowest because we have as- 
sumed Ji < 0. The (0, 0) solution obviously propagates 
as the ferromagnetic state. The solutions {tt, 0), (0, 7r),i.e. 
{I, t, t)i (tj t) 4) pl us their degenerate reversed spin coun- 
terparts can easily be seen to propagate in the up-up- 
down-down stated. 
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FIG. 2: (color online). The classical ground state magnetic 
phase diagram. 



induces the chiral phase which remains stable even in 
presence of a uniform magnetic field^. 

Third case: Now we consider the ID frustrated ferro- 
magnetic model with DM interaction. The solutions for 
general equations arc 

{9,6') = (0,0) where 7 = 
(9,6') = (0,tt) = (tt,0) where 7 = 
(6,9') = (9 ,-9q) where 
1 7 



4o 



cos t) sin f/o 



(11) 



These results show, the ferromagnetic and up-up-down- 
down orders exist only in the absence of the DM interac- 
tion. As soon as the DM interaction increases from zero, 
the ground state of the system goes to the chiral phase, 
independent of the frustrated parameter a (Figf2|). Thus, 
using the cluster LK method, the classical ground state 
phase diagram of the ID spin- 1/2 frustrated ferromag- 
netic model with added uniform DM interaction consists 
of a single phase: "chiral". 



Wc listed the energies for the various stationary solu- 
tions: 

hferro = h c (0, 0) = ^ (-2 + 2a), 
o 

huudd = h c (0, 7r) = -^-(-2a), 

h chira i = h c (9 , -9 Q ) = ^(-J- - 2a). (9) 

8 4a 

By equating these energies in pairs we found only one 
critical point, a c = 0.25. The ground state is in the 
ferromagnetic phase for a < 0.25 and is in the chiral 
phase for a > 0.25. 

Second case: We have put J2 = (or a = 0) which 
means we have just ID isotropic ferromagnetic Heisen- 
berg model plus a uniform DM interaction. In this case 
there is only one solution: 

(9,9') = (9q,—9q) (chiral type) where 

tan 6>o = -7. (10) 

Which means that immediately after turning on the DM 
interaction, the long range chiral order is created in the 
classical ground state phase diagram. In principle, as 
soon as the DM interaction is applied, the ground state 
of the system undergoes a phase transition from the fer- 
romagnetic phase into a chiral phase. 

It is better to emphasis that the induced effects 
of quantum fluctuations on the ID antiferromagnetic 
Heisenberg model with DM interaction had been studied 
theoretically and experimentally^— . From quantum 
point of view, it has been found that the DM interaction 



III. QUANTUM PHASE DIAGRAM 

In this section, to explore the nature of the spectrum 
and the quantum phase transition, we used the Lanczos 
method to diagonalize numerically chains with length up 
to N = 24. 

First, in different subspaces we have computed the 
ground state energy of chains with J\ = — 1 and dif- 
ferent values of the DM interaction as a function of the 
parameter a. The setting of boundary conditions and 
observation of energy with respect to change of param- 
eters is important for precise analysis of wave stated, 
but because of our limitation and size effects, we have 
just set the periodic boundary condition. In Fig. [3J we 
present results of these calculations for the value of DM 
vector D = 0.05 and chain size N = 20. It can be seen 
that the ground state energy for weak DM interaction 
is nearly degenerate up to the first critical frustration 
a ci = 0.24 ± 0.01, which value is obtained by extrapola- 
tion technique. Also degeneracy in the mentioned region 
shows that in the absence of the frustration, a — 0, the 
spectrum of the model is gapless. By turning on the frus- 
tration, spectrum remains gapless up to the first critical 
value of the frustration, a Cl = 0.24±0.01. As soon as the 
frustration increases from a Cl , the ground state is non- 
degenerate and exists in the subspace with total S z = 0. 
Fig. [3] is also good analysis, especially for the proof of 
ferromagnetic sate. 

From the viewpoint of the symmetry our model is 
completely different than the models studied in the 
literature2Sr^a^. The latter models have the Z 2 x 
U(l) symmetry. In this case, the absence of the chiral 
LRO is natural, and the spontaneous breaking of the Z2 
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FIG. 3: (color online). Ground state energy of the system 
in different subspaces versus the frustration parameter a for 
chain with length N = 20 and DM vector D = 0.05. 

chiral symmetry, which resulted in the chiral LRO, is in- 
teresting and was indeed the main topic in those studies. 
On the other hand, in the present model, the Z-i chi- 
ral symmetry is broken by the DM interaction explicitly. 
Therefore, the appearance of the chiral LRO should be 
natural; if the LRO is absent, it should be regarded as 
a surprising result and there must be a very interesting 
mechanism which recovers the symmetry behind the re- 
sult. 

In order to recognize the different quantum phases, we 
have implemented the Lanczos algorithm on finite size 
chains to calculate the lowest eigenstate. The first in- 
sight into the nature of different phases can be obtained 
studying the chiral order parameter. It has been shown 
that the DM interaction can create the chiral phase in 
the ground state phase diagram of the spin system a 32 ' 33 , 
which is characterized by nonzero value of the chirality 

x" = ^ (12) 

3 

where Xj = (Sj x Sj + i) 2 . One should note that there 
are two different types of the chiral ordered phases, 
gapped and gaplcs o 45 ' 47 . In Fig. Ufa), we have pre- 
sented calculated results on the chiral order parameter 
as a function of the frustration parameter a for a value 
of DM vector D = 0.05, including different chain lengths 
N = 16, 20, 24. As is clearly seen from this figure, there is 
no long-range chiral order along the z axis in the regions 
a < a ci and a > a C2 . However, in the intermediate 
region, a ci < a < a C2 , the ground state shows a pro- 
found chiral order in the z direction. Since the study 




n 



FIG. 4: (color online), (a) The chiral order parameter \ z 
as a function of the frustration parameter a for the value of 
DM vector D — 0.05, including different chain lengths N = 
16, 20, 24. In the inset the chiral correlation function is plotted 
as a function of the chain length N for the value of frustration 
a = 0.36. The diverging behavior of the correlation function 
shows that the chiral ordering in the intermediate region is 
true long range, (b) The correlation (xj Xj+n) as a function 
of n for chain length TV = 20 and D = 0.05. 

of correlation functions can give us deep insight into the 
characteristics of the ground state, we define the chiral 
correlation function as 

N 

C~~ =Y.^UUn)- (13) 

n=l 
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It is important that in a phase with true long-range order 
the correlation functions should diverge as N — > oo. To 
check the existence of the chiral long-range order in the 
thermodynamic limit N — > oo of the system, we have 
plotted in the inset of Fig. QJa) the N dependence of C z 
for different values of frustration in the regions a > a Cl . 
As is seen from this figure, only in the intermediate region 
a Cl < a < a C2 , there is a diverging behavior which shows 
that the chiral order in the intermediate region is true 
long-range order. On the other hand, the constant value 
in the region a > a C2 shows that the C z /N takes zero 
value in the thermodynamic limit N — > oo. For grater 
emphasis, we have plotted in Fig. HJb), the correlation 
function (%| Xj+n) as a function of n for a chain length 
N = 20. Clearly is seen that the long distance correlation 
of the chirality only exist in the region a Cl < a < a C2 . 
The mirror symmetry is the result of the periodic bound- 
ary condition. Therefore, by studying the ground state 
energy and the chiral order parameter, we showed that in 
the intermediate region a Cl < a < a C2 , the ID spin-1/2 
frustrated ferromagnetic model with added DM interac- 
tion is in the gapped chiral order phase. 

In what follows we will present our numerical study 
about the presence of dimcr phase. In order to under- 
stand the nature of the dimer phase let us assume we 
have a ID frustrated antiferromagnetic model, where the 
sign of J\ will be changed through the 7r rotations around 
the z axis of the spins on every second sites and the 
model will be transformed to the ID frustrated ferro- 
magnetic model. From the fact that the ground state 
at the Majumdar-Ghosh point (J2 = Ji/2) is given by 
the product of singlet dimmers (^75 (| tl) — I 4-t}))^- 
Through the above 7r-rotation transformation, the dimcr 
unit for J2 = — Ji/2 is replaced by the triplet state 
72 (I tl) + I 4-t}}^- To this reason, to find additional in- 
sight into the different phases we are focused on existing 
of the dimer phase in our ID frustrated ferromagnetic 
model. The order parameter characterizing the dimer 
phase is given as£ 



(14) 



This order parameter shows alternating signs along the 
spin chain. Our results for the DM vector value D = 0.05 
and different chain sizes N = 16, 20, 24 are presented in 
Fig. EJ As can clearly be seen from this figure, in the 
Luttinger liquid region, a < a Cl , there is not any long- 
range dimer order. As soon as the frustration increases 
from the first critical value, the dimer ordering increases 
from zero very rapidly. The oscillations of the dimcr 
order in the intermediate region at finite N are the result 
of level crossing between ground state and excited states 
of the model. On the other hand, overlapping of the 
numerical results in the region a > a C2 , shows a divergent 
behavior of the correlation function of the dimer order 
parameter by increasing the size of chain N. This justifies 
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FIG. 5: (color online). The dimer order parameter, d, as 
a function of the frustration parameter a for the value of 
DM vector D — 0.05, including different chain lengths N = 
16,20,24. 
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FIG. 6: (color online). The quantum ground state magnetic 
phase diagram. 



that the true long-range dimer order exists in the region 
a > a C2 of the ground state phase diagram. 

We have to mention that for other values of the DM 
interaction (D < 0.3), we did the same numerical experi- 
ment and found the same quantum picture of the ground 
state magnetic phase diagram contains of: (I.) gaplcss 
Luttinger liquid phase in the region a < a Cl (II.) gapped 
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chiral phase in the intermediate region a Cl < a < a C2 
and (III.) gapped dimer phase for a > a C2 (FiglHl) - 

IV. CONCURRENCE 

In this section we study the entanglement of forma- 
tion as a measure of the entanglement. Entanglement 
stands some similarity to classical correlation, but it dif- 
fers in some important respects, including the fact that 
entangled objects can violate Bells inequality^. May be 
one of the most characteristic differences is this: if two 
similar quantum objects are completely entangled with 
each other, then neither of them can be at all entangled 
with any other object, but there is no such restriction on 
classical correlations^. 

We compute the entanglement between two sites which 
is known as the concurrenc e) 50 ' 52 

C lm = 2max{o,cQ,C^}, (15) 

where 

eg = vW, + fi) 2 + (95-fi} 2 

where M z is magnetization along z axis and g^ = 
(S^S^) is the correlation function between spins on sites 
/ and m. 

The numerical Lanczos results on the concurrence are 
shown in Fig. [7] From our numerical results we found 
that the concurrence of two spins that are nearest neigh- 
bor is equal to zero in all regions of the quantum ground 
state phase diagram. In this figure the concurrence of 
one spin and its next nearest neighbor is plotted as a 
function of the frustration parameter a for chain lengths 
N = 16, 20, 24 and DM vector D = 0.05. In the Luttinger 
liquid and dimer phases, the correlation between differ- 
ent components, gf^, gf^ is equal to zero and only in the 
chiral phase has non-zero value. As it can be seen from 
Fig. [3 in the Luttinger liquid region, a < a Cl , since the 
correlations gf£, gf^ are equal, therefore the concurrence 
is equal to zero. On the other hand in the intermediate 
chiral region, a Cl < a < a C2 , due to non-zero value of 
the correlations between different components the con- 
currence increases by increasing the frustration and at 
the second critical value a = a C2 is maximum. In sector 
of the dimer ground state, a > a C2 , the correlations gf£, 
dim are ec l ua l an d close to the maximum value ^21. In 
this case, as it is clearly seen from the Fig. the concur- 
rence has the almost constant value that obeys from the 
equation C; m =| gf^+gf^ \— 0.5. The deviation from the 




FIG. 7: (color online). The entanglement between next near- 
est neighbors Cj j+2 as a function of the frustration parame- 
ter a for the value of DM vector D = 0.05, including different 
chain lengths N = 16, 20, 24. The inset shows the concurrence 
as a function of the inverse length 1 /N for different values of 
the frustration a = 0.60,0.72,0.80,0.92,1.00. Finite linear 
extrapolated results show that in the thermodynamic limit 
N — Y oo, the next nearest neighbors are entangled in the 
dimer phase. 



saturation value 0.5 is a result of the induced quantum 
fluctuations by DM interaction. In the inset of Fig. [JJ the 
iV-dependence of the concurrence is investigated for dif- 
ferent values of the frustration in the dimer sector of the 
ground state phase diagram. Finite linear extrapolated 
results show that in the thermodynamic limit N — > oo, 
the next nearest neighbors are really entangled in the 
dimer phase. Thus, in the ID frustrated ferromagnetic 
spin-1/2 model in presence of the DM interaction, only 
the next nearest neighbors are entangled in the gapped 
chiral and dimer phases. 

It is expected that quantum nature of entangle- 
ment might provide explanatory and predictive power 
for the investigation of quantum phase transition s] 51 ' 52 . 
Thus concurrence -pairwise entanglement- which refers 
to quantum correlations has emerged as one of the im- 
portant tools. In our model Eq.([2|), only the next nearest 
neighbors quantum correlation lead to concurrence have 
different nonzero value in the gapped chiral and dimer 
phases and zero value in the Luttinger liquid phase. 

On the other hand, entanglement generation and dis- 
tribution are problems of central importance in perform- 
ing quantum-information(QI) tasks, like teleportation 5 ^ 
and quantum cryptography 5 ^. From the QI perspective, 
it would be attracting to create sizable entanglement be- 
tween particles that arc located at a distance larger than 
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FIG. 8: (color online). The entanglement between end spins 
as a function of a, with D — 0.05 for different chain lengths up 
to TV = 20, with open boundary condition. The inset shows 
Cin as a function of the inverse length 1/7V for different values 
of the frustration a = 0.60, 0.76, 1.00, 2.00. The extrapolated 
results in the thermodynamic limit TV — > oo are less than 
one. 



a few sites. This fact naturally leads to the concepts of 
long distance entanglement (LDE) as a sort of quantum 
order parameter—. There are some models same as the 
dimerizcd- frustrated model, spin-1 Hcisenbcrg chain with 
biquadratic interaction and so on, which able to produce 
LDE££. Here we are interested to check this feature in 
our model. As we showed in previous sections our spin 
frustration model leads to dimerized phase if the param- 
eter a exceeds a certain critical value (a C2 ), so the pro- 
cess of dimerization with the open boundary condition 
enables frustrated model to produce LDE. In Fig. [8j we 
give the numerical Lanczos results on the long distance 
entanglement of frustrated Hciscnberg chain in the pres- 
ence of the DM interaction, particularly D ~ 0.05. As it 
can be seen from Fig.[8j in contrast next-nearest neighbor 
entanglement, the system has got nonzero LDE which in- 
creases rapidly by increasing a up to a C2 and then after 
that it reaches its saturation value about C\-n — 0.7. 
In the inset of Fig. [8j the TV-dependence of the LDE 
is investigated for different values of the frustration in 
the dimer sector of the ground state phase diagram. Lin- 
ear extrapolated results show that in the thermodynamic 
limit TV — > oo, the end-spins are entangled in the dimer 



phase for open boundary conditions. 

V. CONCLUSION 

To summarize, we studied the effect of a uni- 
form Dzyaloshinskii-Moriya (DM) interaction on the 
ground state phase diagram of the one-dimensional (ID) 
isotropic frustrated ferromagnetic spin-1/2 model using 
the analytical cluster method and numerical Lanczos 
technique. Our classical analysis results based on LK 
method show that immediately after turning on the DM 
interaction, the ground state of the system goes to the 
chiral phase, independent of the frustration parameter 
a = -p^j. Thus, we concluded that the classical ground 
state phase diagram of the ID frustrated ferromagnetic 
spin-1/2 model with added uniform DM interaction con- 
sists of a single phase: "chiral". 

To find the quantum corrections of adding the DM 
interaction, we did a very accurate numerical experi- 
ment. We implemented the Lanczos algorithm to find 
the ground state in finite chains. Based on the numeri- 
cal results of the order parameters and correlation func- 
tions, we identified three different phases. In the region 
a < a Cl , the ground state of the system is in the Lut- 
tinger liquid phase. In the region a > a C2 the true long- 
range dimer order exists in the ground state phase dia- 
gram and in the intermediate region a Cl < a < a C2 , the 
system finds in the gapped chiral phase. 

On the other hand, we tried to answer this question: 
"Are the spins entangled in different sectors of the ground 
state phase diagram?" Our numerical results for the en- 
tanglement showed that the nearest neighbor spins are 
not entangled in all regions of the quantum ground state 
phase diagram. But, the next nearest neighbor spins are 
entangled in the gapped chiral and dimer phases, and be- 
came saturate in the region a > a C2 , which can be take 
it as another witness the presence of dimer phase. We 
have even checked this model with open boundary condi- 
tions as an candidate for LDE model, where calculations 
showed that the model reaches its saturated value around 
a C2 which is more than next nearest neighbor entangle- 



ment case Cin > C 
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